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Abstract. In rotor-router aggregation on the square lattice H? ^ par- 
ticles starting at the origin perform deterministic analogues of random 
walks until reaching an unoccupied site. The limiting shape of the clus- 
ter of occupied sites is a disk. We consider a small change to the routing 
mechanism for sites on the x- and y-axes, resulting in a limiting shape 
which is a diamond instead of a disk. We show that for a certain choice 
of initial rotors, the occupied cluster grows as a perfect diamond. 



1. Introduction 

Recently there has been considerable interest in low-discrepancy deter- 
ministic analogues of random processes. An example is rotor-router walk 
|PDDK96] ■ a deterministic analogue of random walk. Based at every vertex 
of the square grid Tl? is a rotor pointing to one of the four neighboring ver- 
tices. A chip starts at the origin and moves in discrete time steps according 
to the following rule. At each time step, the rotor based at the location of 
the chip turns clockwise 90 degrees, and the chip then moves to the neighbor 
to which that rotor points. 

Holroyd and Propp |HP09j show that rotor-router walk captures the mean 
behavior of random walk in a variety of respects: stationary measure, hitting 
probabilities and hitting times. Cooper and Spencer [CS06| study rotor- 
router walks in which n chips starting at arbitrary even vertices each take 
a fixed number t of steps, showing that the final locations of the chips 
approximate the distribution of a random walk run for t steps to within 
constant error independent of n and t. Rotor-router walk and other low- 
discrepancy deterministic processes have algorithmic applications in areas 
such as broadcasting information in networks |DFS08| and iterative load- 
balancing [FGSlOj. The common theme running through these results is 
that the deterministic process captures some aspect of the mean behavior 
of the random process, but with significantly smaller fluctuations than the 
random process. 
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Rotor-router aggregation is a growth model defined by repeatedly releas- 
ing chips from the origin o G Z^, each of which performs a rotor-router walk 
until reaching an unoccupied site. Formally, we set Aq = {o} and recursively 
define 

Am+l = AmU {Zm} (1) 

for m > 0, where Zm is the endpoint of a rotor-router walk started at the 
origin in and stopped on exiting Am- We do not reset the rotors when a 
new chip is released. 

It was shown in |LP081 ILP09j that for any initial rotor configuration, the 
asymptotic shape of the set Am is a Euclidean disk. It is in some sense 
remarkable that a growth model defined on the square grid, and without 
any reference to the Euclidean norm = {xf + x^)^''^, nevertheless has a 
circular limiting shape. Here we investigate the dependence of this shape 
on changes to the rotor-router mechanism. 

The layered square lattice is the directed multigraph obtained from the 
usual square grid by reflecting all directed edges on the x- and y-axes 
that point to a vertex closer to the origin. For example, for each positive 
integer n, the edge from (n, 0) to (n — 1,0) is reflected so that it points 
from (n, 0) to (n -|- 1,0). Only edges on the x- and y-axes are affected. 
Rotor-router walk on is equivalent to rotor-router walk on with one 
modification: if the chip is on one of the axes, and the rotor points along 
the axis towards the origin after it is turned, then the chip ignores the rotor 
and moves in the opposite direction instead. 

For n > 0, let 

Dn = {{x,y) G ^2 : |x| + \y\ < n} . 
We call Dn the diamond of radius n. Our main result is the following. 

Theorem 1. There is a rotor configuration po, such that rotor-router ag- 
gregation {Am)m>o on TI? with rotors initially configured as po satisfies 

^2n(n+i) = Dn for all n > 0. 

A formal definition of rotor-router walk on and an explicit description 
of the rotor configuration po are given below. 

Let us remark on two features of Theorem [T| First, note that the rotor 
mechanism on is identical to that on except for sites on the x- and y- 
axes. Nevertheless, changing the mechanism on the axes completely changes 
the limiting shape, transforming it from a disk into a diamond. Second, not 
only is the aggregate close to a diamond, it is exactly equal to a diamond 
whenever it has the appropriate size (Figure [T]). 

In [KLlOj . we studied the analogous stochastic growth model, known as 
internal DLA, defined by the growth rule ([T]) using random walk on Z^. 
This random walk has a uniform layering property: at any fixed time, its 
distribution is a mixture of uniform distributions on the diamond layers 

Lm = {ix,y) G : |x| -\-\y\=m} , m>l. 
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Figure 1. The rotor-router aggregate of 5101 chips in the 
layered square lattice W' is a perfect diamond of radius 50. 
The colors encode the directions of the final rotors at the 
occupied vertices: red = north, blue = east, gray = south 
and black = west. 

It is for this reason that we call W' the layered square lattice. 

As a consequence of the uniform layering property, internal DLA on 
also grows as a diamond, but with random fluctuations at the boundary. 
Theorem [T] thus represents an extreme of discrepancy reduction: passing to 
the deterministic analogue removes all of the fluctuations from the random 
process, leaving only the mean behavior. For a similar "no discrepancy" 
result when the underlying graph is a regular tree instead of Z^, see |LL09j . 

To formally define rotor-router walk on T?, write ei = (1,0), 62 = (0, 1) 
and let i? = ( ? "j^^ ) be clockwise rotation by 90 degrees. For each site 
z G ^2 \ {0} there is a unique choice of a number j G {0,1,2,3} and a 
point w in the quadrant 

Q = {{x,y) G ^2 : rc > 0, y > 0} 

such that z = Ww. Given j and w = we associate to z = Ww a 

4-tuple {e^, el, el, el) of directed outgoing edges, where 

_ \{z, z + We2) \ii = 2 and x = 0; 
^ \{z,z + R^^^e2) otherwise. 

Thus, for z G Q (hence j = and w = z) the edges 6^,6^,6^,6^ point 
respectively north, east, north, west when z is on the y-axis; and north, east, 
south, west when z is off the y-axis. For z in another quadrant, the directions 
of 6^,6^,6^,6^ are obtained using rotational symmetry. To the origin we 
associate the 4-tuple (6^^, e\, 6^, 6^) where 6^ = (o, i?*e2) for i = 0, 1, 2, 3. 
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Figure 2. Left: The layered square lattice Tl? . Each di- 
rected edge is represented by an arrow; multiple edges on the 
X- and y-axes are represented by double arrows. The origin o 
is in the center. Right: The initial rotor configuration po- 

For every z G Z^, let E^, be the multiset {e^, e;^, e^, e^}. If e = e!, G -E^, we 
denote by e"*" the next element e*"*"^ "^"^ of Ez under the cyclic shift. The 
layered square lattice is the directed multigraph with vertex set V = 'E'^ 
and edge multiset E = E^, where edges that appear twice in Ez have 

multiplicity two (Figure^ left). Thus every vertex has out-degree four, and 
every vertex except for the origin and its neighbors has in-degree four. 

The initial rotor configuration po appearing in Theorem [T] is given by 

po{z) =el zG (3) 

It has every rotor in the quadrant Q pointing north, and rotor directions in 
the other quadrants given by rotational symmetry (Figure [2] right). 

We may now describe rotor-router walk on W' as follows. Given a rotor 
configuration p with a chip at vertex z, a single step of the walk consists of 
changing the rotor p(z) to /^(z)"^, and moving the chip to the vertex pointed 
to by the new rotor p(-z)^. This yields a new rotor configuration and a 
new chip location. Note that if the walk visits z infinitely many times, then 
it visits all out- neighbors of z infinitely many times, and hence visits every 
vertex of W' (except for o) infinitely many times. It follows that rotor-router 
walk exits any finite subset of in a finite number of steps; in particular, 
rotor-router aggregation terminates in finitely many steps. 

2. Strong Abelian Property 

In this section we prove a "strong abelian property" of the rotor-router 
model. Theorem [2| which holds on any finite directed multigraph. To prove 
Theorem [T| we will apply the results of this section to the induced sub- 
graph Dn oil? . 

Let G = {y^E) be a finite directed multigraph (it may have loops and 
multiple edges). Each edge e E -E is directed from its source vertex s(e) to 
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its target vertex t(e). For a vertex v £ V, write 

= {e£ E : s(e) = v} 

for the multiset of edges emanating from v. The outdegree of v is the 
cardinahty of Ey. 

Fix a nonempty subset S C V of vertices called sinks. Let V = V \ S, 
and for each vertex u G y , fix a numbering e^, . . . , e^"~-^ of the edges in E,^. 
If e = el £ Ey, we denote by e+ the next element e*+^™°'^'''' of Ey under 
the cyclic shift. 

A rotor configuration on G is a function 

p:V' ^E 

such that p{v) G Ey for all v eV . A chip configuration on G is a function 

Note that we do not require cr > 0. If a{v) = m > 0, we say there are 
m chips at vertex v; if a{v) = —m < 0, we say there is a hole of depth m at 
vertex v. 

Fix a vertex v £ V . Given a rotor configuration p and a chip configura- 
tion (T, the operation Fy of firing v yields a new pair 



where 



and 



Fy{p,a) = {p,a') 
p'{w) : 



p{w)'^ ii w = v; 
p{w) ii w ^ v; 



a{w) — 1 if w = v; 
a'{w) = < a{w) + 1 if w = t{p{v)~^); 
a{w) otherwise. 

In words, Fy first rotates the rotor at v, then sends a single chip from v 
along the new rotor p{v)~^. We do not require a{v) > in order to fire v. 
Thus if (t{v) = 0, i.e., no chips are present at v, then firing v will create a 
hole of depth 1 at v; if cr(v) < 0, so that there is already a hole at v, then 
firing V will increase the depth of the hole by 1. 

Observe that the firing operators commute: FyFy^ = FyjFy for all v,w e 
V'. Denote by N the nonnegative integers. Given a function 

u:V' ^'N 

we write 

= n ^v^""^ 

where the product denotes composition. By commutativity, the order of the 
composition is immaterial. 
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A rotor configuration p is acyclic if the spanning subgrapli {V, 
no directed cycles or, equivalently, if for every nonempty subset A C V 
there is a vertex v £ A such that t{p{v)) ^ A. 

In the following theorem and lemmas, for functions /, g defined on a set 
of vertices A gV, we write "/ = g on A" to mean that f(v) = g{v) for all 
V £ A, and "/ < 5 on ^" to mean that f{v) < g{v) for all v £ A. 

Theorem 2 (Strong Abelian Property). Let p be a rotor configuration and a 
a chip configuration on G. Given two functions ui,U2 : V — )• N, write 

F^^{p,a) = {p,,ai), i = l,2. 

If o'l = o"2 on V' , and both pi and p2 are acyclic, then ui = U2- 

Note that the equality ui = U2 implies that p\ = p2-, and that ai = 02 
on all of V . For a similar idea with an algorithmic application, see |FL10[ 
Theorem 1]. 

In a typical application of Theorem [2j we take cJi = (T2 = on 1/', and u\ 
to be the usual rotor-router odometer function 

u\{y) = #{1 <j<k : Vj = v} 

where f 1, f2, . . . , ffc is a complete legal firing sequence for the initial configu- 
ration (p, a); that is, a sequence of vertices which, when fired in order, causes 
all chips to be routed to the sinks without ever creating any holes. Provided 
ui{v) > for all v G V , the resulting rotor configuration pi is acyclic: in- 
deed, for any nonempty subset A C V , the rotor at the last vertex of A to 
fire points to a vertex not in A. 

The usual abelian property of rotor- router walk |DF9H Theorem 4.1] 
says that any two complete legal firing sequences have the same odometer 
function. The Strong Abelian Property allows us to drop the hypothesis of 
legality: any two complete firing sequences whose final rotor configurations 
are acyclic have the same odometer function, even if one or both of these 
firing sequences temporarily creates holes. 

In our application to rotor-router aggregation on the layered square lat- 
tice, we take V = Dn and S = L^- We will take a to be the chip configura- 
tion consisting of 2n(n -|- 1) -|- 1 chips at the origin, and p to be the initial 
rotor configuration pQ. Letting the chips at the origin in turn perform rotor- 
router walk until finding an unoccupied site defines a legal firing sequence 
(although not a complete one, since not all chips reach the sinks). In the 
next section, we give an explicit formula for the corresponding odometer 
function, and use Theorem [2] to prove its correctness. The proof of The- 
orem [1] is completed by showing that each nonzero vertex in £)„ receives 
exactly one more chip from its neighbors than the number of times it fires. 

To prove Theorem [2] we start with the following lemma. 

Lemma 3. Let u : V' ^ ¥\ , and write 

F-{p,a) = {p^,a,). 
If a = ai, and pi is acyclic, then u = 0. 
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Proof. Let A = {v £ V : u{v) > 0}, and suppose that A is nonempty. Since 
pi is acyclic, there is a vertex v £ A whose rotor pi{v) points to a vertex 
not in A. The final time v is fired, it sends a chip along this rotor; thus, 
at least one chip exits A. Since the vertices not in A do not fire, no chips 
enter A, hence 

^o-l(f) < ^cr{v) 

contradicting a = ai. Therefore, A is empty. □ 

Theorem [2] follows immediately from the next lemma. 
Lemma 4. Let ui,U2 : V ^ N, and write 

F''^{p,a) = {p^,ai), i = l,2. 
If pi is acyclic and o"2 < cri on V , then ui < U2 on V' . 
Proof. Let 

Then (pi,(Ti) is obtained from (/5, cj) by firing only vertices in the set A = 
{v e V : ui{v) > U2{v)}, so 

(7 < cJi on A'^. 

Likewise, (/>2, C2) is obtained from (p, a) by firing only vertices in A^, so 

<5" < o"2 < o"! on A. 

Thus CJ < (Ji on V. Since J2veV ^(^) ~ X^tieV '^i(^) follows that a = ai. 
Taking 

u = ui — min(ui, ii2) 
in Lemma [sj since F^{p, a) = {pi, ai) we conclude that -u = 0. □ 

3. Proof of Theorem [T] 

Consider again the rotor-router model on the layered square lattice Z^. 
We will work with the induced subgraph D„ of Z^, taking the sites in the 
outermost layer L„ as sinks. 

Recall our notation 

Q = {{x,y) G : X > 0,y > 0} 

for the first quadrant of Z^. We have = {0} U (Ui=o^^Q)' ""^^ere R = 
{i~q) is clockwise rotation by 90 degrees. Fix n, and for z = (x, y) G Dn 
write 

£z = n — \x\ — \y\. 

Consider the sets 

C2 = {ix,y) GQnDn^i : X > 0, y > 2, i^^^y) = 2 mod 4} 
C3 = {{x,y) £QnDn-i : X > 0, y > 1, i^^^y) = 3 mod 4} 
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Figure 3. The rotor configurations p2, Ps, • • • , /07 in the first 
quadrant. The lower left corner is the origin in each picture. 
On the axes, the black arrows correspond to the directed 
edge in ([2]), and open- headed arrows to e^. 



and 



Define n„ : Dn-i — )• N by 



1 



c 



where 



u„{z) 



(4) 



(5) 



2n(n + 1) if z = o; 
+ if 2^0; 

and Iciz) is the indicator function which is 1 for z £ C and for z ^ C. 

Let po be the initial rotor configuration ([s]), and define the rotor config- 
uration pn on Dn-i and chip configuration cj„ on Z)„ by setting 

F""(po,(2n2 + 2n+l)5o) = (p„,cJ„). 

From the formula Q it is easy to obtain an explicit description of p„, and 
to verify that these rotor configurations are acyclic for all n > 1 . Figure [3] 
depicts the rotor configurations p2, ps, ■ ■ ■ , pr in the first quadrant. 



Lemma 5. For all n> 1, we have (Tr. 



1 



Proof. The origin has no incoming edges in W' , so it receives no chips from 
its neighbors. Since Un{o) = 2n? + 2n, the origin is left with exactly one chip 
after firing. The sink vertices Ln do not fire and only receive chips. Since 
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Un{z) = 2 for all z G Ln-i, it follows from Q and (|3| that exactly one chip 
is sent to each sink vertex. 

It remains to show that crn{z) = 1 for each vertex z G Dn-i \ {o}, i.e., 
that the number of chips sent to z by its neighbors is one more than the 
number of times z is fired (that is, 1 + Un{z)). To show this, write 

F<{po, {2v? + 2n + l)5o) = a'^) 

where u'^ is given by ([5|. We will argue that (t'^{z) = 1 and that (Tn{z) = 
(y'n{z). By symmetry, it suffices to consider points z = {x,y) in -Dn-i H Q. 
We argue separately in the two cases x = and x > (on the axis and off 
the axis). 

Case 1: x = 0. Under F"^", the site z fires izi^z + 1) times. If y > 1, its 
neighbor z — e2 fires {iz + l)(^z + 2) times, and from ^ and ([S]) we see that 
it sends a chip to z every even time it is fired. Since {iz + 1)(^2 + 2) is even, 
it follows that z — e2 sends ^{iz + 1)(^2 + 2) chips to z. The same is true 
if y = 1, since then £z = n — 1, and the origin o = z — e2 sends ^n(n + 1) 
chips to z. 

The only other vertices that send chips to z under are its left and 
right neighbors z it ei. Since iz±ei = iz — ^, these neighbors fire £z{iz — 1) 
times. We claim that together they send ^£z{iz — 1) chips to z. To see this, 
note that if we fire these two vertices in parallel, they send one chip to z 
every two times we fire. We therefore conclude that 

a^z) = l{iz + mz + 2) + ^4(4 - 1) - 4(4 + 1) = 1. 

To show that an{z) = (y'j^{z), note first that neither z nor z — e2 is in C 
because x = 0. The right neighbor z + ei might be in C, but since 4(^z — 1) 
is even, the last chip sent from z + ei by F"" does not move to z. The 
left neighbor 2: — ei is in C only if 4-ei =^2 — 1 = 3 mod 4, which implies 
4(4 — 1) = mod 4. Hence if z — ei is in C, the last chip sent from z — gi 
by F"^'^ moves west. It follows that F"" and F"" fire z the same number of 
times and send the same number of chips to z, hence (T„(z) = (t^(z) = 1. 

Case 2: x > 0. To argue that cr^(z) = 1, as an initial step we itnfire 
every vertex on the positive x-axis B = {(m,0) G : m > 0} once. Since 
all initial rotors on B point east, this turns all these rotors north without 
affecting the number of chips at z (nor at any other vertex of Q). 

Now we apply F^". By firing the four neighbors of z in parallel, it is 
easy to see from ([2]) that they send one chip to z every firing round, since 
after every round exactly one of their rotors points to z. Hence, firing 
these neighbors 4(^2 — 1) times each sends ^2(4 — 1) chips to z. Since 
4-ei = 4-e2 = ^2 + 1) the two neighbors z — ei and 2 — 62 each fire 

(4 + i)(4 + 2) - 4(4 - 1) = 44 + 2 

additional times under F"^'-^ . Considering what happens when they are fired 
in parallel shows that they send one chip to z every two times they fire, 
meaning that 21 z + 1 additional chips are sent to z. 
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Finally, to obtain a'^ we must fire every vertex in B once more. But since 
i?^n fires each vertex in B an even number of times, their rotors are now 
pointing either north or south, so firing them once more does not affect the 
number of chips at z. Hence 

a'^iz) = 4(4 - 1) + (24 + 1) - 4(4 + 1) = 1. 

To finish the proof, we now argue that cr„(z) = (t'^{z). First note that 
since 4(4 + 1) is even for all v £ -Dn-ij it follows from Q that the last chips 
sent from z it ei by F"" do not move to z. However, consider the neighbor 
z + e2. If 4+e2 = ^2 ~ 1 = 3 mod 4, its final rotor points north after firing 
4(^2 — 1) times, while if 4 — 1 = 2 mod 4, its final rotor points south. It 
therefore follows from the definition of C, that F"" sends one fewer chip 
from z + e2 to z than F^" in case 4 = 3 mod 4 and y + 1 > 2. Likewise, 
F"" sends one fewer chip from z — e2 to z than F"" in case = 2 mod 4 
and y — 1 > 1. But these are precisely the two cases when z S C, hence -F"" 
also fires z once fewer than F"". Therefore, an{z) = ^"^(-z) = 1. □ 

We remark that the rotor configuration pn is obtained from by cycle- 
popping: that is, for each directed cycle of rotors in p^, unfire each vertex 
in the cycle once. Popping a cycle causes each vertex in the cycle to send 
one chip to the previous vertex, so there is no net movement of chips. Let 

be the acyclic rotor configuration obtained from cycle-popping, and let 

where c„(z) is the number of times z is unfired during cycle-popping. Then 

F<{po,{2n^ + 2n + l)5o) = [pllD^). 

By Lemma [5j we have 

F""(po, (2n2 + 2n + l)5o) = (pn, Id J- 

By the Strong Abelian Property (Theorem [2]), it follows that li^ = Un- In 
particular, Cn = Ic and p'^ = pn- 

Proof of Theorem^ For m > 0, let be the smallest integer such that 
2rm{rm + 1) > rn,. Consider a modified rotor-router aggregation defined by 
the growth rule 

where Zm is the endpoint of a rotor-router walk in and stopped on exiting 
Am n Dr^-i. Define Un : Dn-i — )> N by setting 

Un{z) = # of times z fires during the formation of ^2n(n+i)- 

We will induct on n to show that n„ = n„ for all n > 1. Since Un = Un 
implies A2n(n+i) = ^2n(n+i) = by Lemma [sj this proves the theorem. 

The base case of the induction is immediate: ui = ui = A5o- For n > 2, 
in the induced subgraph D„ of with sink vertices L„ we have 

F""(po, (2n2 + 2n + 1)5„) = (p,, Id J 
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by Lemma [5j On the other hand, 

F^"(po,(2n2 + 2n + l)5o) = 

for some rotor configuration p„ on -D„-i and chip configuration (T„ on D„. 
By the inductive hypothesis, ^2n{n-i) = Dn-i, from which it fohows that in 
the formation of ^2n(n-i-i) from A2n(ra-i); rotor-router walks are stopped 
on exiting Dn-i- Together these facts imply that a„ = 1 on -Dn-i- Moreover, 
since po is acyclic, pn is acyclic (each rotor points in the direction a chip 
last exited). The Strong Abelian Property (Theorem [2]) now gives ti„ = Un, 
which completes the inductive step. □ 
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